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This is a limited review and update of the status of Heavy Flavor Physics. After we review the flavor problem 
we discuss a number of topics: recent puzzles in purely leptonic D and B decays and their possible resolutions, 
mixing in neutral B and D mesons, the determination of IVcfcl and \Vub\ from semileptonic decays, and we 
conclude with radiative B decays. 



1. Introduction: The Flavor Problem 

The Standard Model (SM) of electroweak interac- 
\ tions^ correctly accounts for all known particle physics 
data. Hints of small anomalies exits but none is firmly 
' established. If there is new physics (NP) it must be 
hiding, and one good way to hide it is by making it ac- 
tive only at shorter distances than we have yet probed. 
Future high energy particle collision experiments may 
directly probe such new short distance physics. We 
can ask, in the mean time, what are the indirect ef- 
fects of such new physics at the longer distances that 
are probed in current experiments? A model indepen- 
dent way to address this question is by supplementing 
the Lagrangian of the SM with local terms, or "op- 
erators," of dimension greater than four. Such terms 
render the theory non-renormalizable. Hence a mo- 
mentum cut-off Anp must be introduced and the the- 
ory is considered only as an Effective Field Theory 
(EFT) valid only at energies below Anp. Colliding 
particles with center of mass energy in excess of Anp 
surely produces new states that require further spe- 
cific modification of the theory. 

The scale Anp also serves to make the terms in the 
i Lagrangian dimensionally correct. An operator of di- 
' mension n > 4 appears with coefficient c/AJ^p*, with c 
a dimensionless constant. The natural expectation is 
that c is of order unity. A very large number renders 
the theory in-effective, breaking down at energies be- 
low Anp. This does not happen: one simply chooses a 
smaller number for Anp. On the other hand, the co- 
efficients of some operators could be very small. But 
short of explaining why some coefficients unexpectedly 
small, we must assume that in fact we underestimated 
Anp. Hence, we may proceed by assuming c ~ 1 and 
see what current data implies for Anp . 

The EFT generically contains AF = 2 FCNCs, that 
is, terms that induce neutral interactions that change 
flavor by two units. For example one may include 

— ^ [ci((iL7''SL)(rfL7M^i) + C2(wL7^CL)(ui7^CL)] . 

Anp 

(1) 



^Supplemented with right handed neutrinos and masses. 



Then, if one ignores the SM contribution, neutral me- 
son mixing data gives [ij 

^(data) ^ ^ Q Q33^^ ^ (^-^^ ' , (2) 

cr^U(5.9 + 1.0^)xlO-^(^)^ (3) 

There being no reason to expect a cancellation be- 
tween the SM and NP contributions, the NP contri- 
butions should be smaller than c\'^'^^^^ . Therefore Anp 
should be larger than the electroweak scale by some 
four orders of magnitude! 

This in itself is not a problem. But there is one 
good reason to expect NP at the electroweak scale. 
In the SM there are quadratically divergent radiative 
corrections to the higgs mass. In terms of our cut- 
off EFT, the shift in the higgs mass from L-loops is of 
order Anp / (47r)^ so a counterterm must be fine tuned 
to one part per mil to cancel this at one loop, and 
further fine tuned to one per cent at two loops, etc. 
If instead Anp ~ 1 TeV there is no need for any fine 
tuning. This is the Flavor Problem, that NP at the 
EW scale requires extraordinarily small dimensionless 
couplings Ci. 

Much of the work on Heavy flavor physics aims at 
testing the SM in the flavor sector with high preci- 
sion. It gives additional restrictions on NP, that can 
be described as bounds on additional coefficients of 
higher dimension operators. The purpose of this talk 
is to present some of the main results in heavy flavor 
theory. While it is interesting to investigate models of 
NP that address the Flavor Problem, a prerequisite is 
to understand the restrictions that heavy flavors place 
on the models. In the absence of glaring anomalies in 
the data, this is best done by verifying consistency of 
the SM to as high precision as possible. I will focus on 
precision SM determinations, but will here and there 
indicate implications on models of NP. 
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2. Purely Leptonic Decays 
2.1. The evanescent fo, puzzle 

The theory of purely leptonic decays is simple, 

T{Ds tvi) = ^/l, G|m2|KsP(l - ml/ml f 
an 

(4) 

with fo^ the Ds decay constant and V the Kobayashi- 
Maskawa matrix. There are analogous formulas with 
obvious modifications when replacing B, Bg or D for 
Dg. Last year a discrepancy became apparent in the 
value of fo, obtained from Monet-Carlo simulations 
of QCD on the lattice and the one from the experimen- 
tal measurement of the purely leptonic branching frac- 
tion. From a recent compilation of lattice results 

fo = 206(4) MeV, fo, = 243(3) MeV (5) 

while experimentally[l, [J 

fo = 205.8(8.5)(2.5) MeV, fo = 275(16)(12) MeV 

(6) 

While this anomaly was not firmly established, the 
agreement between lattice and experiment in the value 
of f]j suggests the discrepancy in f^^ may well remain 
once the errors are reduced. 

Perhaps for this reason several groups have looked 
for a viable interpretation of this result in terms of 
NP[i,|l,0|- In the SM this is a Cabibbo allowed, tree 
level decay. Hence for the NP to have a significant ef- 
fect it must be neither loop nor Cabibbo suppressed. 
Moreover the mass M of the new particle mediating 
this interaction cannot be too large: for constructive 
interference the amplitude should be about 6% of the 
SM's, so roughly M « M/VOM = 320 GeV. Do- 
brescu and Kleinfeld argued that (i) s-channel charged 
higgs exchange could explain the effect, with <ti jjc 
and both jjc and yr of order unity, but then found 
this explanation disfavoured by D decay data (ii) t- 
channel charge -1-2/3 leptoquark exchange could also 
account for the data but is disfavored by the bound 
on T — !■ ssfi (iii) w-channel charge —1/3 leptoquark 
exchange (d-squark like object, d) is a viable explana- 
tion. They introduce the interaction 

-Clq = Ki{cLil - SLviL)d + n'iCRe^d + h.c., (7) 

which is already present in supersymmetric exten- 
sions of the SM without R parity, and show that for 
l^'g/Kil <C mimc/m^ the resulting interference is au- 
tomatically constructive. In addition, if \ki\ « \kt-\ 
the deviations in fii^ and tv are approximately equal. 

Earlier this year the CLEO collaboration published 
new results on Dg purely leptonic decays to both ^i/ 
and TP final states [8|. Their value for the decays con- 
stant, fo, = 259.5(6.6) (3.1) MeV has significantly re- 
duced errors but also has moved significantly in the 
direction of eliminating the anomaly. 
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Figure 1: Tension in the CKM global fit from B —y tv^. 
The cross corresponds to the experimental values with 1 
sigma uncertainties. The shaded region is the result of the 
global fit performed without using these measurements. 

2.2. B^Tv 

The direct determination of the branching fraction 
for B ^ TV and sin 2/3, which measures the CP asym- 
metry in interference between mixing and decay in 
B J/tpKs, is in slight disagreement with a global 
fit to these quantities from other measurements. This 
is presented in Fig. [T] which shows both the result of 
the direct measurements (as a cross) and that of the 
global fit (as a shaded area) . 

That there is a correlation between these quantities 
can be seen from the following formula, obtained from 
the expression for the B'^ — B'^ mass difference^ Am^, 
by eliminating the decay constant fs in favor of the 
purely leptonic branching fraction: 

BijB ^ Tv) _ - 1 /sin^y 

Awd ~ ^'m\S{xt)\VudV Bb^ \sm-fj 

(8) 

Here S{xt) is an Inami-Lim function with the top 
quark mass as argument, and Bbj is the "bag" con- 
stant that parametrizes the matrix element of the 
short distance four-quark operator that gives rise to 
mixing. The main uncertainties come from Bb^ and 
the angles a and 7. 

The deviation is 2.4 sigmas if one compares the in- 
direct fit prediction for the B tv branching frac- 
tion with the direct measurement. Alternatively one 
can use the fit to determine Bb^ and compare with 
the value determined form Monte-Carlo simulations 



^See next section. 
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Figure 2: The shaded area gives the predicted ratio of 
B semileptonic branching fractions into r vs. electron, 
as a function of r = tan/3/Mj^+, in a two higgs doublet 
model [13]. The dashed-dotted lines give experimental la 
bounds. 

of QCD on the lattice. Then the deviation is 2.7 sig- 
mas. 

As in the Dg decay case this is a tree level process, 
but in contrast, this is Cabibbo suppressed, which 
gives a bit more leeway in giving a NP interpreta- 
tion. In a two higgs doublet model tree level charged 
higgs exchange can easily modify the leptonic branch- 
ing fraction[l3, [HI- In this model the Cabibbo sup- 
pression is still present but there is an amplification 
factor proportional to tan^ /? — {v2/vi)'^, where vi 
and V2 are the vacuum expectation values of the higgs 
doublets responsible for up and down masses, respec- 
tively. The parameter that controls the NP correction 
is r = tan/3/M/f+. The semileptonic branching frac- 
tion will also be affected, and a sensitive test of this in- 
terpretation is the ratio of semileptonic decay branch- 
ing fractions into t versus electrons [T^: see Fig. O 



3. Neutral Meson Mixing 
3.1. Generalities 

The recent measurements of mixing of neutral D 
mesons and their unfamiliar properties suggest we 
begin our discussion by a taking a general look at 
neutral meson mixing. Let us briefly review the 
quantities that enter the description of neutral me- 
son mixing. Parameters p and q are introduced to 
express the physical states in terms of flavor eigen- 
states: \Pl,h) = p\P^) ± q\P^)- These, together with 



the mass and width differences. Am = mn — 
and AF = Th ^ F^, are given in terms of the off- 
diagonal elements of the 2x2 non-hermitian hamilto- 
nian M - |F by 

(Am)2-i(AF)2 = 4|Mi2p-|Fi2p, 
AmAF = 4Re(Mi2Fj2), 
f_ ^ 2Mi*2 - iTl^ 
p2 2Mi2 - iTi2 ■ 

We expect M12 rather than F12 to be more prone to 
modifications from new physics because F is largely 
given by long distance physics. Finally, the decay am- 
plitudes are 

Af^{f\n\P') and Af = {f\n\P'). (10) 

CP violation (CPV) can be searched through pro- 
cesses that probe different quantities. If 7^ 1 
there is CPV in decays, while \q/p\ 7^ 1 gives CPV 
in mixing. A non- vanishing imaginary part of A/ = 
{q/p){Af /Af) gives CPV in the interference between 
mixing and decay. The phase ^12 = arg(— M12/F12) 
is sensitive to NP that may show up in M12. The 
parameter that controls the di-lepton asymmetry is 
ImFi2/Mi2 = (1 - \q/p\^)/{l + \q/p\^); it is non- 
perturbative and hence difficult to compute (the OPE 
is no better that for lifetimes, perhaps worse). 

The behavior of the mixing system depends rather 
sensitively on which of Am and AF is largest. Con- 
sider first the case Am ^ AF which is the situa- 
tion for B and mesons. This condition corre- 
sponds to small Fi2/Afi2 and one can find approx- 
imate solutions of (III): Am = 2|Mi2|(l -I- • • • ) and 
AF = — 2Fi2 cos 012 (1 + •••); where the ellipsis indi- 
cate corrections of order F12/M12. Keeping in mind 
that (j)i2 is suppressed in the SM, we see that the ef- 
fects of NP can only reduce the magnitude of AF. On 
the other hand since q/p = — arg(Mi2)(l + • • • ) time 
dependent CP asymmetries are sensitive to NP that 
may show up in M12. 

The other extreme case has AF ^ Am and this 
condition corresponds to small A/12/F12. As before, 
one can find approximate solutions of Am = 

2|Mi2C0S(/)i2|(l •••), AF = t2Fi2(1 + ---) and 
q/p = — arg(Fi2)(l-l-- • • ) depends weakly on M12 (the 
ellipsis now indicate corrections of order M12/F12). 
For example, if D-mesons satisfy AF 3> Am and there 
in neghgible CPV in the decay, then a.igXx+K- oc 
|Afi2/Fi2p sin(20i2). Hence there is reduced sensitiv- 
ity to NP in M12, even for dominant NP. 

3.2. B°B° and B3s 

CPV in B decays and mixing is discussed at length 
in other talks at this conference. In order to avoid un- 
necessary duplication we limit ourselves to discussing 
the mass difference measurements and theory. 
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In the SM and mixing have the same un- 
derlying dynamics (double W exchange, with virtual 
top-antitop quark intermediate state). Hence the ex- 
pressions'^ for Ato are virtually identical, except for 
obvious change of parameters. Taking the ratio not 
just simplifies the expressions but also cancels some 
uncertainties. Solving for the ratio of KM elements 



td\ 



I Ams 



where ^ 



Bb, f%^ 



(11) 

Here fs and Bb are the decay and bag constants (see 
([8])) introduced in previous sections. The quantity ^, 
required to extract a value for the ratio of KM ele- 
ments, contains all of the hard to estimate hadronic 
physics and has the property that it is unity in the 
SU{2)y symmetry limit (to^ — rns)- Monte-Carlo 
simulations of lattice QCD give 



C = 1.205(52) 
C = 1.258(25)(21) 



FNAL/MILC]]^ 
HPQCD[li] 



(12) 



Using Belle and BaBar accurate measurements of 
Amd and CDF and DO measurements of Am^, Evans 
reports 



\Vtd\ 

\Vts\ 



0.2060 ± 0.0012(exp)+°:°°^j^(theor). 



I hasten to point out that Evans uses a value for ^ 
reported in 2003[16J, one of the early unquenched cal- 
culations. Luckily, this value interpolates between the 
two more recent results, ^ — 1.210(^35), and the er- 
rors are also similar. Rather than fussing over the best 
central value and theoretical error let us stop to think 
how well we can trust lattice calculations with this ex- 
traordinary, ~ 3% precision. In fact, the precision of 
the calculation is only about 16% since only the devia- 
tion of from unity needs be computed. Some other 
rather crude methods should therefore work rather 
well too. For example, a very early computation of 
the chiral logs gives ^ — 1.14[l7]. If we estimate the 
errors by A^^ ^ (mx/A^)^ « 20%, we would write 
^ = 1.14±0.08, or a 7% error. Remarkably, this crude 
determination agrees with the lattice result, within ex- 
pected errors! 



3.3. D°D° 

In the two years that followed the first evidence for 
D"D" oscillation by BaBar 18] and Belledl, these 
two collaborations have reported a number of related 
and improved measurements^!^ i2flfc [2l|, [l^, [13, [2^ 
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Figure 3: Contours of allowed values of the parameters 
X = Am/r and y — AF/F from an HFAG analysis of the 
Belle and BaBar data on D" mixing and decays [27|. 



l26j . It is customary to introduce x = Am/F and 
y = AF/F. The experimental situation is summarized 
in Fig. [3l which shows the allowed region in the x, y 



^Which are rather involved, witness Eq. ISl l. 



The two- IF exchange graph that gives AF = 2 pro- 
cesses in the SM has, in the DD case, an intermediate 
qq' state with q, q' = d, s, b. Since these are light com- 
pared to the W, GIM suppression is very effective in x 
and y. Therefore x, y are small. But they are not per- 
turbatively calculable. They are probably dominated 
by the same long distance physics so x may well be 
comparable in size to y. Moreover, y has a Cabibbo 
suppression factor sin^ 9c and vanishes in the SU (2)„ 
limit (m^. = md). 

The experimental measurement x ~ y ~ 1% is very 
compatible with SM expectations. However, since we 
can't compute accurately, precision tests of the SM are 
not possible. Still one can use the measured values of 
X and y to constrain models of NP. M12 is more sensi- 
tive to NP because, in the SM, it starts at 1-loop while 
F12 starts at tree level. Moreover, NP is short distance 
dominated. Excluding a possible perverse cancellation 
between the short and long distance contributions to 
M12 we can then restrict the NP by demanding that 
its contribution be no larger than what is measured. 
We have already alluded to this in a model indepen- 
dent way in ([3]). The implications for many specific 
models have been studied, including the MSSM and 
a sequential fourth generation of quarks [28|, |2§| . The 
MSSM can evade the bound by the same mechanisms 
that are available for the other flavor problems of that 
model: either make the SUSY breaking scale uncom- 
fortably high or find a mechanism, like gauge medi- 
ation, that implements Minimal Flavor Violation in 
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Figure 4: Results of the t combining all available moment 
data, from Ref. |4Q]. The Ax = 1 contours are shown for 
the t with and without B —> Xs'j data. The analysis is 
performed assuming a IS (bottom) subtraction scheme. 

the MSSM. Because it must be very heavy, a fourth 
generation quark in the two-W exchange graph does 
not have an effective GIM mechanism. But its contri- 
bution can be suppressed by small KM (4 x 4) matrix 
elements, \Vub'Vcb' \ 0.001 for nib' - few lOO's GeV. 



4. Determination of |Vcb 



The method of moments gives a very accurate de- 
termination of \Vcb\ from inclusive semileptonic B de- 
cays. In QCD, the rate dT{B Xciv)/Ax<ly = 
|K6p/(a;,y), where x and y are the invariant lepton 
pair mass and energy in units of ms, is given in terms 
of four parameters: |V"cfc|, a^, TOc and mf,. \Vcb\, which 
is what we are after, drops out of normalized mo- 
ments, (z") — J dx dyf{x,y)z" / J dxdyf[x,y) where 
z — X ov y. Since as is well known, the idea is to 
fix rric and rrib from normalized moments and then 
use them to compute the normalization, hence deter- 
mining \Vcb\- In reality we cannot solve QCD to give 
the moments in terms of rric and m^, but we can use 
the combined HQET/OPE to write the moments in 
terms of nic, rrib and a few constants that parametrize 
our ignorance [il H [H, H 113, HI. These con- 
stants are in fact matrix elements of operators in the 
HQET/OPE. If terms of order 1/toq are retained in 
the expansion one needs to introduce five such con- 
stants; and an additional two are determined by me- 
son masses. All five constants and two quark masses 
can be over-determined from a few normalized mo- 
ments that are functions of E^cut, the lowest limit of 
the lepton energy integration. The error in the de- 
termination of |V"cb| is a remarkably small 2%[33; see 
Fig.m But even most remarkable is that this estimate 
for the error is truly believable. It is obtained by as- 
signing the last term retained in the expansion to the 
error, as opposed to the less conservative guessing of 
the next order not kept in the expansion. Since there 
is also a perturbative expansion, the assigned error is 
the combination of the last term kept in all expan- 
sions, of order [3oal, Q;sAQCD/"ifc and (Aqcd/'716)^. 



This section and the next are concerned with the 
precision determination of the magnitudes of the two 
KM elements Vcb and Vub ■ It goes without saying that 
a precise determination is fundamental in testing the 
KM theory of CP violation in the SM. Traditionally 
these magnitudes are determined from the semilep- 
tonic decay rate of B mesons. Hence these sections 
are largely about the theory of semileptonic B decay, 
and we separate each one into inclusive and exclusive 
sub-sections since the theory is different. 

4.1. Inclusive 

The theory of inclusive semileptonic B decay is 
based on a double expansion, combining the heavy 
mass expansion of Heavy Quark Effective The- 
ory (HQET) with the Operator Product Expansion 
(OPE)0, [Sll, HI. The decay rate can then be ex- 
pressed as an expansion, F = Fq -I- Fi -|- • • • , where 
F„ ~ (AqcD/wb)"- The main results are that Fq is 
actually given by the perturbative expression for the 
semileptonic decay rate for an unconfined h quark and 
that Fi = 0. 



There is only one assumption in the calculation that 
is not fully justified from first principles. The mo- 
ment integrals can be computed perturbatively (in 
the XjrfiQ expansion) only because the integral can 
be turned into a contour over complex energy E away 
from the physical region[33|- However, the contour 
is pinned at the minimal energy, -Ecut, on the real 
axis, right on the physical cut. So there is a small 
region of integration where quark-hadron duality can- 
not be justified and has to be invoked. Parametrically 
this region of integration is small, a fraction of order 
A/mg of the total. But this is a disaster because this 
is parametrically much larger than the claimed error 
of order (A/mg)'^. However, this is believed not to 
be a problem. For one thing, the fits to moments as 
functions of i?cut are extremely good: the system is 
over-constrained and these internal checks work. And 
for another, it has been shown 'Zl'l that duality works 
exactly in the Shifman-Voloshin (small velocity) limit, 
to order 1/toq. It seems unlikely that the violation 
to local quark-hadron duality mainly changes the nor- 
malization and has mild dependence on i?cut , and that 
this effect only shows up away from the SV limit. 
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4.2. Exclusive 

The exclusive determination of \Vcb\ is in pretty 
good shape theoretically, and only last year has be- 
come competitive with the inclusive one. So it 
provides a sanity check, but not an improvement. 
The semileptonic rates into either D or D* arc 
parametrized by functions J-"*, of the rapidity of 
the charmed meson in the B rest-frame, w. Luke's 
theorem [4^ states = JF, = 1 + ©(Aqcd/^c)^ at 
w = 1. The rate is measured at w > 1 and extrapo- 
lated to w = 1. The extrapolation is made with a first 
principles calculation to avoid introducing extraneous 
errors [431 . The resulting determination of \Vcb\ has 
a 4% error equally shared by theory and experiment; 
the theory error is dominated by the uncertainty in 
the determination of J^, J^, at w = 1. FNAL/MILC 
combines the 2008 PDG average for |Vc6|.F'*(l) with 
their computed value = 0.921(13)(20) to ob- 

tain, after applying a small electromagnetic correc- 
tion, \Vcb\ = (38.7±0.9exp±1.0thcor) X 10-3(11. 

There is some tension between theory and experi- 
ment in these exclusive decays that needs attention. 
The ratios of form factors Ri^2 are at variance from 
theory by three and two sigma respectively [i^. Also, 
in the heavy quark limit the slopes of J- and J-^, 
should be equal. One can estimate symmetry viola- 
tions and obtains [46| — ~ 0.19, while experi- 
mentally this is —0.22 ±0.20, a deviation in the oppo- 
site direction. This is a good place for the lattice to 
make post-dictions at the few percent error level that 
may lend it some credibility in other areas where it is 
needed to determine a fundamental parameter. 



5. Determination of |Vub| 

The magnitude \Vub\ determines the rate for B — > 
Xuii'- The well known experimental difficulty is that 
since \ Vub\ ^ \Vcb\ the semileptonic decay rate is dom- 
inated by charmed final states. To measure a signal 
it is necessary to either look at exclusive final states 
or suppress charm kinematically. The interpretation 
of the measurement requires, in the exclusive case, 
knowledge of hadronic matrix elements parametrized 
in terms of form-factors, and for inclusive decays, un- 
derstanding of the effect of the kinematic cuts on the 
the perturbative expansion and quark-hadron duality. 

5.1. Inclusive 

This has been the method of choice until recently, 
since it was thought that the perturbative calculation 
was reliable and systematic and hence could be made 
sufficiently accurate. However it has become increas- 
ingly clear of late that the calculation cannot be made 
arbitrarily precise. The method uses effective field 



theories to expand the amplitude systematically in in- 
verse powers of a large energy, either the heavy mass 
or the energy of the up-quark (or equivalently, of the 
hadronic final state). One shows that in the restricted 
kinematic region needed for experiment (to enhance 
the up-signal to charm-background) the inclusive am- 
plitude is governed by a non-perturbative "shape func- 
tion," which is, however, universal: it also determines 
other processes, like the radiative B Xg'j- So the 
strategy has been to eliminate this unknown, non- 
perturbative function from the rates for semileptonic 
and radiative decays. 

Surprisingly, most analysis do not eliminate the 
shape function dependence between the two processes. 
Instead, practitioners commonly use parametrized fits 
that unavoidably introduce uncontrolled errors. It is 
not surprising that errors quoted in the determination 
of \Vub\ are smaller if by a parametrized fit than by 
the elimination method of [47|. The problem is that 
parametrized fits introduce systematic errors that are 
unaccounted for. 

Parametrized fits aside, there is an intrinsic prob- 
lem with the method. Universality is violated by sub- 
leading terms [ill in the large energy expansion ("sub- 
leading shape functions" ) . One can estimate this un- 
controlled correction to be of order asA/rrib, where 
A is hadronic scale that characterizes the sub-leading 
effects (in the effective theory language: matrix el- 
ements of higher dimension operators). We can try 
to estimate these effects using models of sub-leading 
shape functions but then one introduces uncontrolled 
errors into the determination. At best one should use 
models to estimate the errors. I think it is fair, al- 
beit unpopular, to say that this method is limited to 
a precision of about 15%: since there are about 10 
sub- leading shape functions, I estimate the precision 
as \/TOasA/mb. This is much larger than the error 
commonly quoted in the determination of |Kif,|. 

This is just as well, since the value of \Vub\ from 
inclusives is in disagreement not only with the value 
from exclusives but also with the global unitarity tri- 
angle fit. You can quantify this if you like, but it is 
graphically obvious from Fig. [5l The location of the 
apex of the unitarity triangle differs in the two pan- 
els, and the agreement would be much better if the 
green ring, whose radius is given by iKib/Kbl and is 
dominated in the fit by the determination from the 
inclusive decay, had smaller radius. 

5.2. Exclusives 

The branching fraction Br{B irii/) is known (49j 
to 8%. A comparable determination of \Vub\ requires 
knowledge of the B ^ tt form factor f+{q^) to 4%. 
There are some things we do know about /+: (i)The 
shape is constrained by dispersion relations [50j . This 
means that if we know /+ at a few well spaced points 



Proceedings of the DPF-2009 Conference, Detroit, MI, July 27-31, 2009 



7 



= 1 : 




rl sin2p 










a _z 


= , ,« 


/ ^ \ 






Figure 5: Latest CKMfitter fit of data on the p — rj plane describing the unitarity triangle|9|]. The left panel includes only 
measurements of CPV angles, while the right panel excludes them. The green ring, |Vuij|, is dominates by the inclusive 
semileptonic determination. 



we can pretty much determine the whole function /+ 
(ii)We can get a rough measurement of the form fac- 
tor at = rn^ from the rate for B 7r7r[5l|. Thij 
requires a sophisticated effective theory (SCET) anal- 
ysis which both shows that the leading order contains 
a term with /+(m^) and systematically characterizes 
the corrections to the lowest order SCET. It is safe 
to assume that this determination of /+(mj) will no1 
improve beyond the 10% mark. 

Lattice QCD can determine the form factor, at leas1 
over a limited region of large q^. The experimental 
and lattice measurements can be combined using con- 
straints from dispersion relations and unitarity[53]. 
Because these constraints follow from fundamentals, 
they do not introduce additional uncertainties. They 
improve the determination of \Vub\ significantly. The 
lattice determination is for the g^-region where the 
rate is smallest. This is true even if the form fac- 
tor is largest there, because in that region the rate 
is phase space suppressed. But a rough shape of 
the spectrum is experimentally observed, through a 
binned measurement [49| . and the dispersion relation 
constraints allows one to combine the full experimen- 
tal spectrum with the restricted-q^ lattice measure- 
ment. The best lattice calculations are in good agree- 
ment; however use the same MILC ensembles [52|. 
They give 3.55(25)(50) and 3.38(36) for lO^lKbl- The 
11% error in \Vub\ is dominated by lattice errors. 



5.3. Alternatives 

Exclusive and inclusive determinations of | Vub \ have 
comparable precisions. Neither is very good and the 
prospect for significant improvement is limited. Other 
methods need be explored, if not to improve on exist- 
ing \Vub\ to lend confidence to the result. A lattice- 
free method would be preferable. A third method, 
proposed a while ago[5J|, uses the idea of double 
ratios |55] to reduce hadronic uncertainties. Two inde- 
pendent approximate symmetries protect double ra- 
tios from deviations from unity, which are therefore 




Figure 6: CKMfitter constraints in the p~ rj plane includ- 
ing the most recent inputs in the global CKM fdSX[. The 
\Vub\ constraint has been split in two contributions: from 
inclusive and exclusive semileptonic decays (plain dark 
green) and from B'^ — > tv (hashed green). The red hashed 
region of the global combination corresponds to 68% CL. 



of the order of the product of two small symmetry 
breaking parameters. For example, the double ratio 

(/s.//bJ/(/d,//dJ = {fBjfoJ/ifBjfD,) = 1 + 

0{ms/mc) because fsJlBa = foj foa = 1 by SU {'d) 
flavor, while fBjfo, = fsJfDa = ^/mjrnb by 
heavy flavor symmetry. One can extract \Vub/VtsVtb\ 
by measuring the ratio. 



dTjBq ^ piiy)/dq^ 
dr{Bd^ K*i+i- )/dq^ 



^tsVtb] 



87r2 
■ a2 



iV(g2 



■Rb, 
(14) 

where q^ is the lepton pair invariant mass, and N{q'^) 
is a known function 56]. When expressed as functions 
of the rapidity of the vector meson, y = Ev/my, the 
ratios of helicity amplitudes 



Rb — 



Rd = 



(15) 

are related by a double ratio: RB{y) = RD{y){^ + 
0{ms{m~^ — m^^))). This measurement could be 
done today: CLEO has accurately measured the re- 
quired semileptonic D decavsfSTl [ssj. 
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More methods are available if we are willing to use 
rarer decays. To extract \Vub\ from Br(_B+ —t r^Vr) = 
(1.8 ± 0.6) X 10~''[5^ one needs a lattice determina- 
tion of fs- This was discussed in Sec. 12.21 although 
the determination is still imprecise and relies on lat- 
tice hadronic parameters, it gives an even larger value 
for \Vub\, see Fig. IH] which, however, uses as input 
Br(B+ ^ T+Ur) = (1.73 ± 0.35) x 10-^ with sub- 
stantially lower errors than in the combined BaBar- 
Belle result of Ref. [s^ . Since we want to move away 
from relying on non-perturbative methods (lattice) to 
extract \Vub\ we have proposed a cleaner but more dif- 
ficult measurement, the double ratio 



r(_B„ 



\VtsVtb\^ 



fsl fh 
Id/ fi 



(16) 



In the SM Br(B, -> fi+n') « 3.5 x 10"^ 
x(/B^/210MeV)2(|T43|/0.040)2 is the only presently 
unknown quantity in the double ratio and is expected 
to be well measured at the LHC|60j. 

The ratio r{B+ T+v)/T{Bd fJ-'^fJ-") gives us a 
fifth method. It has basically no hadronic uncertainty, 
since the hadronic factor fs/fBa = Ij by isospin. 
It involves|14&P/|VtdVt{,p, an unusual combination of 
CKMs. In the p — rj plane it forms a circle centered 
at ~ (—0.2, 0) of radius ^ 0.5. Of course, measuring 
T{Bd ^JL^ n^) is extremely hard. 

In a sixth method one studies wrong charm decays 
Bd,s DX (really bq u p) . This can be done 
both in semi-inclusive decays [6l| (an experimentally 
challenging measurement) or in exclusive decaysj62l| 
(where an interesting connection to Bd^s mixing ma- 
trix elements is involved). 



6. Rare Radiative B decays. 

The rare decays B Xg-f and B Xgi^i^ are 
flavor changing neutral processes that occur first at 
one loop in the SM. As such they are sensitive probes 
of NP. They are complimentary to other FCNCs and 
put stringent bounds on NP models beyond what is 
obtained from neutral meson mixing measurements. 
Here I will focus on the radiative decay, B Xgj 
both because of its higher rate and because of recent 
progress in theory. Both the total rate and CP violat- 
ing asymmetries can probe NP. The average of exper- 
imental measurement of the rate is rather precise[63j. 



Br(B 



E^>1.6 GeV 



(355 ± 24+^^0 ± 3) X 10" 



(17) 

Here we have indicated that the measured rate is 
only for energetic photons, > 1.6 GeV. The com- 
bination of data requires some mild extrapolation 
since measurements have differing photon energy cuts. 
BaBar recently reported a rather strong constraint in 



the CPV asymmetry, 0.033 < Acp{B Kj) < 0.028 

The theory B ^ Xg^ has two parts. The first one 
is the computation of the low energy effective Hamil- 
tonian. This is necessary in order to re-sum large 
logarithms, lii(mt/mb) or ln(Mn//mb, in the pertur- 
bative expansion. The second step is the computation 
of the rate from this effective Hamiltonian. 

The effective Hamiltonian, to lowest or- 
der in an expansion in Gp, is Ti = 

-{iGF/V2)V^lV^,J2LlC^i^^)Q^■ The are 
dimension 6 AB = —AS* = 1 operators. Roughly 
they are the tree level four quark operator and the one 
it mixes pronouncedly with, Qi = {s ^CL){cLl^bL) , 
Q2 = {sLltiT°-CL){cLl^T°-bL), four penguin opera- 
tors, Qs-e = ab)^q{cp'al) (with r matrices in 
color and spinor space), and two transition magnetic 
moment operators Q7 = {emb/16n^)sLO''^'^bjiF^^ and 
Qs = {ggmb/l6Tr^)sLcr>"'T''bB.G1^^- The problem is to 
compute reliably the coefficients at a low renormaliza- 
tion scale, fj, ~ rub- This requires computation of the 
coefficients at a short distance scale, /i ~ Mw, and 
then using the renormalization group to "run" the 
coefficients to the physical scale, ^ nib, for which 
the anomalous dimensions of the operators needs to 
be computed. The leadin g lo garithms (LO) were first 
summed 20 years agof65i l66j. The result is a correc- 
tion of more than 30% to the un-resummed coefficient 
Crinib) (and therefore a whopping 60% effect in the 
rate). To achieve accuracy comparable with present 
experimental measurement it is important to re-sum 
the next-to- leading (NLO) and the next-to-NLO 
(NNLO) logs. This is a challenging enterprise that 
has taken the better part of two decades. The 
NNLO calculation requires two loo p m atching of 
Ci„6^67'|, three loop matching of C7-«[68j|. three-loop 
calculation of the (1-6) x (1-6) and (7-8) x (7-8) 
blocks of the anomalous dimension matri x |69| . and 
four-loop of the (1 - 6) x (7 - 8) blockjlj^ The 
magnitude of the NNLO coefficients are, roughly, 
|Ci,2(mb)| ^ 1, |C3,4,5,6(mh)| < 0.07, Cjinib) ~ -0.3 
and Cs{mb) ~ —0.15. 

The second step, the computation of the rate from 
this effective Hamiltonian is no smaller a feat. To 
match the accuracy of the NNLO coefficients one 
needs a perturbative calculation of the matrix ele- 
ments of Qi-fito three-loops [iJl and of Qts to two- 



loops [72, l73|, |7j|. The calculation is not complete, 
some interference terms are missing; see Ref. jTSj 
for a detailed account and for an account of non- 
perturbative effects. 

The SM prediction for the branching fraction re- 
stricted to energetic photons > 1.6 GeV is[76| 

Br(B ^ X,7)|^>i.6 Gov = (3.15±0.23) x lO'^ (18) 

in magnificent agreement with the experimental aver- 
age Eq. 1171). 
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Figure 7: Direct and indirect bounds on A/^± in the Two 
Higgs Doublet Model, type II, as a function of tan /3, from 
Ref . [t^] ■ The shaded areas are excluded by the constraints 
at 95% CL. 



The effectiveness of this process in limiting models 
of NP is nicely illustrated in Fig. [7| which shows the 
constraints form a variety of measurements on a two 
higgs doublet model of type II (charge-2/3 quarks get 
masses from the vacuum expectation value (VEV) of 
one higgs, charge — 1/3 quarks form the VEV of the 
other higgs) in the (Af//± , tan/3) plane. Several mea- 
surement nicely compliment each other, but it is clear 
that radiative B decays plays a leading role in exclud- 
ing parameter space. 

We close with a couple of remarks on B X££. 
The NNLO calculation in the SM requires one less 
loop than B so it has been complete for quite 

some time. Much attention has been given to the 
observation that the forward-backward asymmetry in 
B K*U has a zero in the SM[7i. The presence 
and location of the zero suffer little from hadronic un- 
certainties and from contamination from non-resonant 
B KttU decays [T^. Less well known is the fact 
that at large invariant lepton-pair mass B — > K*£i is 
well understood and that tests of the SM can be done, 
largely free of hadronic uncertainties, by a method of 
double ratiosfi^l as described in Sec. 15.31 Unfortu- 
nately experimentalists have not conducted this test, 
even if the data is available. 
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